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• Concluding Remarks 


Presentation Objective: 

Apply a synergistic multiscale modelling technique to woven composites to 
determine underlying reasons for nonlinear response 

• Understand influence (i.e., primary, secondary, etc.) of architectural parameters 
(e.g., fiber/void volume fraction, weave geometry, tow geometry, void geometry)at 
multiple length scales on the mechanical response of CMCs. 

• Analyze the significance of effects and compare to material scatter 
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Model Efficiency 


Goal is to Balance Efficiency vs Fidelity 
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NASA 's Integrated multiscale Micromechanics 

Analysis Code (ImMAC) Suite 




MAC/GMC is Evolving Anisotropic Thermoelastic 
Inelastic and Damage Constitutive Model 
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Fidelity vs. Efficiency in Composite Micromechanics 

Comparison of Local Stress Invariants 

Transverse Loading; 50% Glass/Epoxy 
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Individual Stress Components 
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Failure Criterion for Strength and Durability 


Subcell Elimination Criterion 
•Max. Stress Theory 
•Max. Strain Theory 
•Tsai-Hill Theory 
•Tsai-Wu 
►SIFT 

•Elastic Allowables 


Progressive Damage Criterion 
• Scalar Damage (Tri axial) 
• MMCDM 

•Smeared Crack Band 
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Integrated Multiscale Analysis of 
Arbitrary Composite Structures with 


FEAMAC 


Synergistic Multiscale Modeling 

• Embed micromechanics within FEAat 
element integration points 

• New tool for micro/macro analysis of 
composite structures: 

FEAMAC 

• Localize/homogenize on the fly 
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Utilize Novel Multiscale Generalized Method of Cells 
(MSGMC) For Concurrent Analysis of 

Wn ve>n/Rr aided Hnmnnsite Systems 
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Problem Definition 

Weave Properties 
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5 Harness Satin 
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Type 

5HS 

Overall Fiber Volume Fraction 

36% 

Tow Volume Fraction 

78% 

Tow Width 

1.25 mm 

Tow Spacing 

0.34 mm 

Thickness 

2.5 mm 

Matrix 

CVI-SiC 
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Tow Properties 
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Fiber Vol Fraction within Tow 
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Multiscale Generalized Method of Cells(MSGMC) 

Overview 


• Newly developed recursive GMC methodology 

• Each length scale in each subcell can call a separate GMC analysis 


• Works for any arbitrary multiphase material 
• Elastic / Inelastic / Damage 




Multiscale analysis can determine 
local stresses at different length 

scales 
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Macroscale (Weave) Two Step Homogenization 

To compensate for lack of normal-shear coupling within GMC a two-step 
homogenization scheme is employed for woven composites. 

(Bednarcyk & Arnold, IJSS, 41 , 2003) 
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Constituent Constitutive Model and Strain Localization 

Microscale 



Fiber/ 
Interface 


Assume Fiber and Interface Linear 
Elastic Hashin Fiber Failure Criteria 
( 1980 ) 
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Assume Linear Elastic with a Scalar 
Damage constitutive relationship 


^{apy){py\ _ ^ j _ ^{ocPr}{Pr} j ^{^Pr}{Pr} £ {aPr}{Pr} 

Matrix damage driven by magnitude of triaxialitv 


°H 

S H 



^ °^H > G critical 

f = 3 £ H nK — cj H — 0 



nAs% x + /L'4 +1 

(a4 +1 + 4 1 ) 

(V +1 = s l + A^ /+1 ) 


K° = initial bulk modulus 



Full Multiscale Modeling of 5HS Weave with 

Porosities 
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5HS and most other orthogonal weaves can be discretized 
into 8 unique subcell groups. Furthermore model tow and 
matrix with voids using lower scale RUCs 


Matrix with 
void 
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Three Void Modeling Schemes Considered 

Voided Matrix Response Achieved via Separate GMC Analysis 

No Voids Evenly Distributed Voids Localized Voids 
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Gold = 12.7% voids ; Red = 90% voids; Blue = 5% voids 
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Localization of porosity significantly influences 
failure response 

a) Knee - 33% delta 

b) Strain to Failure - 25-55% delta 

Assuming uniform distribution of voids similar to no voids 
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Simulation Identifies Local Damage Events / Mechanisms 
Explaining Nonlinearities in Macro Stress Strain Curve 


Assuming 5HS RUC with 



Fiber: Elastic , Hashin Fiber Failure Criteria (includes 
shear stress) 

Interface: Elastic (very compliant 1/20 th ) 

Matrix: Elastic, Hydrostatic-Driven Damage 
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Study Effects Of Micro, Meso, And Macro 
Parameters on Macroscale Response 



Architectural Parameter 

Relevant Length Scale 

Values 

Tow Fiber Volume Fraction 

Meso 

0.46,0.48,0.50 

Tow Void Volume Fraction 

Meso 

0.01,0.05,0.07 

Tow Aspect Ratio 

Macro 

8,10,12 
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Influence of Varying Matrix Material Parameters 

on the Macroscale Response 
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Depicts Entire Range Of Macro Response Curves Given the 
27 Variations In Architectural Parameters 

Utilized Localized Void Model 


Architectural Variations 
clearly contribute to 
variation in measured 
material response. 

• Initial Modulus » 24% 

• UTS *2% 

• 1 st matrix cracking » 
16% 

• Post matrix cracking 
Modulus » 24% 

• 8 f impacted » 16% 




Strain 
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Probability 


Assumed Normal Distributions for Architectural 

Parameters 

Normal Distribution Probability Plot* 
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Note: Material Properties held fixed 
at Baseline Values; Void shape - 
sheet like 


Parameter 

Mean 

Std. Dev. 

Tow Fiber Volume 
Fraction 

0.48 

0.033 

Tow Aspect Ratio 

8 

.533 

Tow Void Volume 
Fraction 

0.05 

0.01 

Localized Weave Void 

0.75 

0.05 


Volume Fraction 


*Bonacuse, P., Subodh M., and Goldberg, R.; ’’CHARACTERIZATION OF THE AS MANUFACTURED VARIABILITY IN A CVI 
SIC/SIC WOVEN COMPOSITE, Proceedings of ASME Turbo Expo 2011, GT2011, June 6-10, 2011, Vancouver, Canada 
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Procedure for Incorporating Stochastics 
Requires Significant Computation Resources 



Weave Type 

Time / Increment 
(sec) 

Typical 

Increments 

Total Time 
(sec) 

No. of 
Subcells 

5HS (1x1) 

12 

200 

4000 (1.1 hrs) 

93,800 

PW (1x1) 

1.5 

150 

225 

18,840 

PW (6x6) 

53 

200 

10600 (2.9 hrs) 

678,276 



Macro Stress-Strain Response Curves Given 
Stochastic Assumption of Architectural Parameters 

Utilized Localized Void Model 



"determined from normality assumption 
using bilinear approximation 


ASH 


Arizona State University 



Secant Through Thickness Moduli (E zz ) Degrades 
With Loading As Does In-plane (E^) 



Note: In composites many material “properties” 
evolve with loading history ! 


E zz ; GPa 


Normal in ^ Normal out 
E zz = 59.02 ± 12.5 ; 68.2% confidence 


E zz = 59.02 ± 25 ; 95.4% confidence 
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Loading Histories with Unloading Are Critical For 
Deducing Mechanisms Driving Nonlinear Response 



Examine Plain Weave Discretization to Study 
Architectural Parameters on Structural Scale 


Tow Aspect Ratio= 
width/thickness 
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Subcell group properties determined 
from lower length scales 



Arizona State University 






Macroscale - Plain Weave Discretization 

Assumes Normal Distribution for all Architectural Parameters 


600 


Macroscale (1x1) 


Statistical Reconstruction 95% Confidence" 


ro 400 



58 Cases Evaluated 


0.2 0.4 0.6 

Strain (%) 


ST 400 



PW slightly less stiff and more nonlinear than 5HS 

rm 


train (%) 

*determined from normality 
assumption using bilinear 
approximation 
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Stress (MPa) Stress (MPa) 


Sensitivity To Architectural Features Changes With 
Increasing Structural Scale: Plain Weave 



4x4 RUC 
8 cases 


0.4 

Strain (%) 


0.8 
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Sensitivity To Architectural Features Changes With 
Increasing Structural Scale: Plain Weave 

Structural Scale (6x6) Statistical Reconstruction 95% Confidence* 




Property 

Ell 

(GPa) 

PLS 

H (GPa) 

ct uts 

Mean 

211 

111 

73.7 

460 

± a 

9.5 

10 

6.5 

42.5 

±2a 

19 

20 

13 

85 
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Comparison of Reconstructed 95% (2a) Confidence 
Plain Weave Stress-Strain Response 

Blue = 1x1, Yellow = 6x6 



■ Composite Stiffness, PLS (first 
matrix cracking), Secondary 
Modulus, statistically unaffected 
by increasing size of RUC 

■ Failure stress/strain is the only 
value that we can say with 95% 
confidence is influenced by 
architectural details 


Property 

Ell (GPa) 

Ell (GPa) 

PLS 

PLS 

H (GPa) 

H(GPa) 

a UTS 

C UTS 

1x1 

6x6 

1x1 

6x6 

1x1 

6x6 

1x1 

6x6 

Mean 

209.5 

211 

116 

111 

74.5 

73.7 

512.5 

460 

± a 

13 

9.5 

10 

10 

6.5 

6.5 

15 

42.5 

± 2a 

26 

19 

20 

20 

13 

13 

30 
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Conclusion 

1. Demonstrated that a synergistic analysis using the multiscale generalized 
method of cells (MSGMC) can accurately represent woven CMC tensile behavior 

(loading/unloading) 

• 4 level of scales analyzed 

• Nonlinear behavior due to damage - demonstrated by unloading 

• Critical invariant is l 1 (brittle) not J 2 (metals) 

• Failure mechanisms capture via local continuum damage model 

2. Non-uniform distribution of voids/porosities must be incorporate within the RUC 
- accurate deformation and failure response 

3. Variations in Weave Parameters (micro, meso, and macro) appear to contribute 
to variation in measured material macrolevel response. 

a) Primary Variables appear to be 

• Constituent material constants (micro) 

• Spatial distribution of void locations (meso); shape is sheet like 

b) Secondary Variables appear to be 

i. Tow void content (meso) 

/'/. Tow Aspect Ratio (meso) 

Hi. Tow volume fraction (macro) 

4. Assuming Normal Probability Distributions — ► showed that only the ultimate 
failure stress/strain (statistically speaking) is influenced at the structural level 
by lower scale features . 
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Future Work 


1. Examine the influence of these parameters on the time-dependent 
material response and corresponding life. 

2. Incorporation of constituent property distribution in the analysis 

3. Incorporate environmental degradation (due to oxidation / moisture) 

4. Multivariate statistics and stochastic processes for coupled 
architectural/material parameters 

5. Incorporate MSGMC into ImMAC 5.0 
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FMC Stress (MPa) Modulus (GPa) 


Average Values of Four Key Composite Response 

Attributes: E, PLS, H and g uts 

Avg. Modulus Avg. Secondary Modulus 




Avg. First Matrix Cracking Stress Avg. Fail Stress 



RUCs RUCs 


FMC Stress (MPa) Modulus (GPa) 


Average Values of Four Key Composite Response 

Attributes: E, PLS, FI and g uts 

Remember 5x5 has lowest DoF 



Std. First Matrix Cracking Stress 
12 , 



RUCs 


Std. Fail Stress 
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